function out = sspa(input)

inputPha = angle(input);
inputMag = abs(input);

p = 2;
p2 = 2 * p;
A0 = 1;

temp = (1 + (inputMag / A0).^(p2));
temp = nthroot(temp, p2);
MagDis = inputMag ./temp;

out = MagDis .* exp(1i * inputPha);
end